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ABSTRACT 

In this paper we investigate the nonlinear development of the most unstable 
Gortler mode within a general three- dimensional boundary layer upon a suitably 
concave surface. The structure of this mode was first identified by Denier, Hall &z 
Seddougui (1991) who demonstrated that the growth rate of this instability is 0(G* ) 
where G is the Gortler number (taken to be large here), which is effectively a measure 
of the curvature of the surface. Previous researches have described the fate of the 
most unstable mode within a two -dimensional boundary layer. Denier Hall (1992) 
discussed the fully nonlinear development of the vortex in this case and showed that 
the nonlinearity causes a breakdown of the flow structure. 

The effect of crossflow and unsteadiness upon an infinitesimal unstable mode was 
elucidated by Bassom &: Hall (1991). They demonstrated that crossflow tends to 
stabilise' the most unstable Gortler mode, and for certain crossflow/frequency combi- 
nations the Gortler mode may be made neutrally stable. These vortex configurations 
naturally lend themselves to a weakly nonlinear stability analysis; work which is de- 
scribed in a previous article by the present authors. Here we extend the ideas of 
Denier and Hall (1992) to the three- dimensional boundary layer problem. It is found 
that the numerical solution of the fully nonlinear equations is best conducted using a 
method which is essentially an adaption of that utilised by Denier and Hall (1992). The 
influence of crossflow and unsteadiness upon the breakdown of the flow is described. 


Research was supported by the National Aeronautics and Space Administration under NASA con- 
tracts Nos. NAS] -18605 and NASI 19480 while the author wa.s in residence at the Institute for Computer 
Applications in Science and Engineering ( ICASE), NASA Langley Research Center. Hampton, VA 28665 
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§1 Introduction 

The aim of this article is to further our understanding of the effects of unsteadiness 
and crossflow upon the fully nonlinear development of unstable Gdrtler modes. The 
initial derivation of the governing equations for these modes was given by Gdrtler 
(1940) whose insults were modified by Hammerlin (1956). These early works ignored 
nonparallel effects present due to boundary layer growth and Smith (1955) added 
terms in an attempt to rectify this deficiency. Until recent years it was unclear as to 
the importance of the nonparallel terms; this question was resolved bv the results of 
Hall (19S2 a,b. 1983). In this series of papers Hall showed, using a combination of 
asymptotic and numerical techniques, that for order one wavenumber vortices there is 
no unique neutral linear stability curve. More precisely, the stability characteristics of 
such wavenumber modes are entirely dependent upon the initial form and location of 
the disturbance. However for small wavelength vortices a unique neutral curve does 
exist and on this curve the vortex wavenumber k is O(G^) wdiere G is the (assumed 
large) Gdrtler number. For further reading concerning the development of the stability 
theory for Gdrtler vortices the reader is referred to Hall (1990). 

The most unstable Gdrtler mode (i.e. that infinitesimally-sized vortex which has 
the largest growth rate) was obtained by Denier, Hall k Seddougui (1991) and Tim- 
oshin (1990). By considering the stability properties of the essentially viscous modes 
of wavenumber (){ G * ) together with those of the inviscid modes of 0(1) wavelengths 
it was possible to identify an intermediate wavenumber regime in which the vortex 
growth rate is largest. These unstable modes are confined to a region of depth of 
0{G~s) and possess growth rates of size 0(G*). The stability properties of 0(0*) 
wavenumber vortices are deduced by solving a sixth-order ordinary differential system 
with appropriate boundary conditions and the solution of this svstem reveals that the 
unique most unstable mode has wavenumber k = 0.476Gs and growth rate 0.312G*. 

To date there has been relatively little attention paid to the nonlinear stability 
properties of Gortler vortices. Perhaps the first work was performed by Aihara (1976) 
who attempted to describe the nonlinear evolution of Gortler vortices using parallel 
arguments. Later calculations by Hall (1988) remedied these defects and showed that 
as 0(1) wavenumber modes evolve downstream so the energy of the flow 7 concentrates 
itself in the fundamental and mean flow correction. This suggests far downstream the 
flow 7 can be adequately described by a mean field/first harmonic structure and such a 
configuration w 7 as elucidated for short wavelength modes using both weakly nonlinear 
and fully nonlinear approaches by Hall (1982 b) and Hall k Lakin (198S) respectively. 

A fuller description of the nonlinear stability of Gdrtler modes may be found in 
Denier Hall (1992). In that paper the authors argued that in a number of practical 
systems, especially where significant curvature occurs such as the case of flow 7 over 
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turbine blades, one would expect the small localized surface imperfections may well 
trigger the most unstable linear Gortler mode. (This conclusion relies on the result 
of the linear receptivity theory given by Denier, Hall k Seddougui (1991).) This 
motivated a careful study of nonlinear evolution of the most unstable mode which was 
tackled in the following way. Denier k Hall (1992) took an arbitrary form of initial 
disturbance at a typical streamwise location, say x = 1. By integrating the linear 
form of the vortex equations over a long distance, up to x — 101. a flow profile was 
obtained which is dominated by the most unstable mode described by Denier et al 
(1991). From .r = 101 onwards the full nonlinear equations replaced the linear ones 
and the most unstable vortex mode was marched further. Typically 8 or 16 harmonics 
of the fundamental were retained during this calculation. It was found that at a critical 
point the flow’ contains a region of reverse flow and the analysis is then no longer valid. 
Denier k Hall (1992) interpreted this breakdown as being responsible for the vortices 
moving away from the wall and into the core of the boundary layer. 

The effect of crossflow and unsteadiness on the most unstable Gortler mode was 
discussed by Bassom k Hall (1991). The primary result arising from this work was 
the demonstration that a relatively small crossflow could completely stabilise the most 
unstable mode. Additionally, by allowing for vortex unsteadiness, it was shown that 
suitable combinations of crossflow, vortex frequency and wavenumber could lead to 
neutrally stable configurations. A weakly nonlinear stability analysis pertaining to 
such configurations was conducted by Bassom k Otto (1992) who derived classical 
‘Stuart-Watson' (1960) type evolution equations for near-neutral modes. They con- 
cluded that the weak-nonlinearity has a stabilising effect and derived equations for 
the supercritical equilibrium amplitudes. 

The results of Denier k Hall (1992) and Bassom k Hall (1991) provide the mo- 
tivation for the current study. Within a two-dimensional boundary layer the effect of 
nonlinearity on the most unstable mode tends to lead to a finite-distance breakdown 
whereas crossflow' appears to stabilise the flow. With these two mechanisms tending 
to have opposite effects it is clear that in many practical situations, where three- 
dimensionality is undoubtedly important, it is of great interest to determine which 
of these two conflicting behaviours dominate. We attempt to answer this question by 
considering the full nonlinear vortex equations and employ numerical techniques which 
are similar to those used in Denier k Hall (1992) but modified in certain ways (detailed 
later). These improvements significantly speed up the computations and allow us to 
obtain a greater range of results than those found by Denier k Hall (1992). 

The structure of the remainder of this article is as follows: in section 2 the fully 
nonlinear equations are derived and a brief description of the numerical procedures 
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used are outlined in section 3. Details of the results are presented in section 4 and in 
section 5 some brief conclusions are drawn. 


§2 Formulation 

As in Hall (1985) we consider a boundary layer flowing over the cylinder i) = 0, 
— oc < z < x where the r-axis is a generator of the cylinder, y measures the distance 
normal to the surface and x denotes the distance along the surface. The Reynolds 
number R e and Gortler number G are defined by 


Re 


UqL 

v 


G = 2 R}8, 


where [ 0 is a typical flow velocity in the x -direction, L is a characteristic streamwise 
lengthscale and u is the kinematic viscosity of the fluid. Furthermore the curvature of 
the cylinder is supposed to be ±\ 0 (f ) and with these definitions 8 = L/b, where b is 
a typical radius of curvature of the cylinder (Gortler vortices are typically observed in 
flows over concave surfaces which corresponds to the choice \o >0 ). The Reynolds 
number is supposed to be large whilst 8 is sufficiently small so that as 8 — ► 0 the 
parameter G is fixed and is of order one. The basic three-dimensional boundary layer 
is taken to be of the form 


u = U 0 (u ( X. Y ),R e *v (X, Y ),R e h*w( X, Y )) (l + O (r7 1 ) ) 

l 

where X = x/L and Y = yR? / L and the crossflow parameter A* is of order one. 

It is convenient to define the scaled spanwise coordinate Z = zR] /L and let T he 
the temporal variable scaled on L/Uq . The basic velocity profile is perturbed by the 
quantity 

(u ( T, X , Y,Z), R7 * V ( T, X, Y.Z). R~ * W ( T. X, Y. Z )) , 

_ 1 

and the pressure field by R e * P(T,X,Y,Z). Substituting this flow form into the 
continuity and Navier-Stokes equations yields the system 


Ux + by + \\ z = 0, 

—Ut 4- Uyy + Uzz — «y C — uU.x—uxU — vUy — \ m wU s 

= UU X + VUy + WU Z , 

■Vt + Vyy + Vzz — G\uU — Py — tfV'y — vx U — vVy — VyV — \*wV : 


(2.1 a) 

( 2 . 1 h) 

(2.1c) 
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— W'f + Wyy T W"zZ ~ Pz ~ u\\ \ — A *W \V — t'W y — A*V ii’y — A U’H ; 

= uwx + vwy + ini'z, 


(2. hi) 


where terms of relative order Q(R, 2 ) have been neglected. It is worth noting at this 
point that the linearised system studied by Bassom k Hall (1991) is obtained by setting 
the right-hand sides of ( 2.1 b-d) to zero whereas the nonlinear equations examined 
by Denier k Hall (1992) to determine the development of steady nonlinear vortices in 
two-dimensional boundary layers can be retrieved by setting A* and Oj equal to zero. 

We now invoke the scalings proposed by Denier e.t. al. (1991) who demonstrated 
that in high Gdrtler number flows the most unstable vortices have 0(G * ) wavenumbers 
and are confined to a layer of thickness of 0(G~$) adjacent to the cylinder. These 
modes have a spatial growth rate of 0{G *) and we use the results of Bassom k Hall 
(1991) who illustrated that the three-dimensionality of the basic flow significantly 
affects the two-dimensional stability results once the scaled crossflow parameter A* 
becomes 0(G%). Therefore it is convenient to define the 0(1) crossflow parameter A 
bv 

A* = G*\. ( 2.2 a) 

To reflect the fact that the vortices are confined to a region of depth 0(G ~ * ) adjacent 
to the cylinder we introduce the 0( 1 ) boundary layer coordinate y defined by 

y = Gh\ ( 2.2 b) 

and in this layer the basic flow may be expanded as a Maclaurin series of the form 

u = G~* fi\\(X)y + 5 i*n{X )y 2 + \G~s fii3(X)y 3 + , (2.2c) 

w = G~* H 2 \(X)y + \G~ T H 22 {X)y 2 + ~G~ r /.t 2 ^{X)y :i + ■•• . (2.2 d) 

To determine the form of the vortex disturbance we appeal to the findings of Bassom 
k Otto (1992) who identified the crucial perturbation size at which the governing 
equations become fully nonlinear (although these authors made no attempt to solve 
these fully nonlinear forms). The disturbance forms are then 

U = G~* fi r o + G~ * u x 4 . G~h T 2 + ■ • ■ j , V = G* (Vo + G 5 Vj + G 5 V 2 + •••'). 

( 2.3a, b) 

W = G 7 (Wq + G~ * W\ + G ~ * W 2 H ) , P = G 7 > (Pb 4- G~ 5 P\ + G 5 P 2 + • • ■) , 

(2.3c, d) 

where U 0 , V’ 0 , H'o, P 0 , U \ , . . . etc. are all functions of X , y, the temporal and the span- 
wise variable. It is now convenient to implement the result of Bassom k Hall (1991) 
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that the leading order behaviour in the downstream coordinate is purely oscillatory 
and we do this by introducing the coordinate and temporal variations given by 


A' = G 


3 

5 


Z = dx + G 


l 

5 - 


T = G“t/. 


2.4</ 


where d = 
become 


-V' 2 i / ii- This then implies that the streamwise and spanwiso denvativ<'s 


d 

dX 


G* 


dx 


- dCA 


d_ 

ol' 


0 

oz 




d 

d~z 


(2Ah) 


The desired governing equations are obtained by substituting (2.2)-(2.4) into equations 
(2.1). Leading order terms in the momentum equations yield that H' 0 = dU 0 . Next 
order terms in equation (2.1 b) give the first relationship 


( d 2 d 2 y 2 a d d d \ r , 

\d^ + OZ 2 ~ ~d~z ~ dT ~ f ' uy fa ) U ° ~ fiuV ° = RHSl - f 2 - 5 ") 

where the precise form of RHS t is given presently and a = A // 22 - d/r 21 . A second 
equation is derived by following a procedure similar to that described by Bassom X Hall 
(1991). By considering the second order terms in the y and r momentum equations 
(2.3), eliminating the pressure by cross-differentiation and applying the continuity 
equation it is a routine but lengthy task to obtain 


d 2 d 2 y 2 a d d 

dy 2 + 0z 2 2~dz ~ dT 


t*n y- 


d 


dx ) \dy 


d 2 


+ 


d 2 

dz 


Vo+o 


dV 0 

dz 


•\o/HU/ 


d 2 U u 

dz 2 


RHS . 


(2.56) 

where again we shall specify RHS 2 shortly. The forms of RHS 2 includes reference to 
the combination \\ i — 3U\ which therefore needs to be expressed in terms of quantities 
with subscript zero. This is best accomplished by integrating the continuity equation 
to give 


W\ 


$U\ = <f>o(x,y) - J 


'dUp dV 0 ' 

dx dy 


dz. 


(2.5c) 


Tlie aim of this article is to consider the nonlinear evolution of modes which are periodic 
in the spanwise direction, with fundamental wavenumber (\oMn ) 5 ^ sa }% and periodic 
in time. Now it is advantageous to introduce scalings first proposed by Bassom k Hall 
(1991). It is found that if we transform according to 


d„ 


(xo^n ) 5 kd y , d : -» (xo/x?i)* kd : , d z -» 


t'o — Xo 5 h n 6 - Vo — \ o I 1 ] i V. 


2 1 - 
<^0 = Xo Hurt' 


A = 


d T -► XoUhdr, 

^ (^ 11^22 — /* 2 l /<12 ) 


2 / f n Xo 5 


(2.6a — c) 
(2 M-f) 


( 2 . 6 ( 7 -/,) 
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then the transformed equations (2.5 a,b) are independent of the partictdar boundary 
layer under consideration. Therefore the ensuing results are potentially relevant to 
a wide class of basic configurations. In this nondimensionalised coordinate system 
we restrict, ourselves to flow quantities the fundamental mode having spanwise and 
temporal variation given by exp [/ (c + (it)], so that (l represents the nondimensional 
frequency of the fundamental vortex component. 

Upon making the transformations (2.6) the leading order vortex equations (2.5 
a.b) become 


r , ° 2 ° 2 

L ( — + 


dy 2 dz 2 


2 A dV y d 2 U 1 <9 2 S (,) 1 d 2 S {2) 1 d 2 S (?) 

v + k*~dz ~ P Hz 2 ~ k 2 ~dz 2 ~Y 2 dzdy ~ Y dxdy ' 




where the operator L is defined according to 


d 2 xb 3 2 0 Ay 2 1 y 

= ajy + aP" ~ ~ Par ~~ n a7 


(2.76) 


(2.7c) 


and 


50. = c ^ + ,v^ + 4 

ax ay a^ 






dU_ dV 
dx + dy J 




+ 7 


r2 


s 


( 2 ) 


4> 


. a 


+ K i 


- r 

k J. [ dx dy 


dU , dV 
+ *- 


dz 


H( 

du r t _ i f r 

fc Jz [ dx 


dU , dV 

!h + k ~d y 


S ^ = U^ + kV^ + k 

dx dy dz 


+ k 


dV 
dy J 


(2.7rf) 


z.,e 


(2-7/) 


To complete the governing equations it is necessary to determine 0. This is achieved 
by noticing that W\ — BU\ satisfies 

dP 

L ( IT ] — 3d i ) = ay \ o H — 1* S^ 2 ^ . 

Identifying the component of this equation independent of z yields 

2n 


2 7T 


d 2 


y 


d 


[dy 2 k*dx\ 


<P 


= f/H 


+ S (2) dz 


h7 g) 
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To dose the system requires specification of appropriate boundary conditions. Clearly 
we require the velocity components £'. T' to vanish on y - 0 (and by continuity so 
must <91 fdy). Additionally the mean flow terms o = 0 on tj — 0 and in order that 
the disturbance' be confined to the boundary layer we demand that the streamwise 
velocity l tends to some function of .r as y — > oc. 

We notice that equations (2.7) are the appropriate generalisations of those' solved 
by Denier anel Hall for the most unstable nonlinear vorte'x within a two elimensional 
boundary layer (their equations are' recovered by setting A = o = 0 and setting <9y=0). 
This allowed us to compare numerical results against their previously published ernes 
as a check of our numerical methods. 

§3 Numerical Methods 

The methods employed to solve system (2.7) are similar in essence to those used 
in Denier Ac Hall (1992). However it was found necessary to introduce' a number of 
amendments to their code in order to speed up the computations which in turn had 
the benefit of allowing an increased number of results to be obtained. 

Denier Hall (1992) solved the two-dimensional counterparts of equations (2.7) 
as follows. They decomposed each of the flow velocities UA'.W into their Fourier 
components and rewrote the governing equations in terms of these components. By 
utilising a scheme based upon that implemented by Hall (19SS) they obtained finite 
difference equations for the mean flow and harmonic terms. These equations contain 
only one streamwise derivative and a straightforward method, based upon solving one 
tridiagonal and one pentadiagonal system, may be used to march the solution from 
one streamwise station to the next. For more details of the practicalities of the scheme 
the reader is referred to Denier Hall (1992). 

One marked difference between our present work and that of Denier <L’ Hall ( 1992) 
is that we chose to calculate the nonlinear terms in physical space rather than trans- 
form space, as this is a computationally ‘cheaper’ policy and so allowed us to retain 
more modes in our calculations. This necessitates transforming from Fourier space 
to physical space, effecting the calculations and extracting the Fourier coefficients. It 
was decided to employ Fast Fourier Transforms to do this which has the benefit of 
reducing the cost of the transformations from 0( N 2 ). (the cost of reducing N modes 
an d using the Fourier transform directly) to 0(5iV logiV — 6iV). The code used was 
based on the original Cooley and Tukey (1965) algorithm and was thus limited to .V 
being an integral power of 2. 

Two other changes were made to the code used by Denier & Hall (1992). As 
was described in the introduction, these authors were concerned with the development 
of the most unstable Gortler mode within a two-dimensional boundary layer. To 



achieve this most unstable mode Denier Sz Hall (1992) took an arbitrary disturbance 
and then marched the linear stability equations for a long distance downstream before 
changing to the nonlinear system. This of course had the effect of ensuring that the 
most unstable component of the original disturbance was then dominant over all the 
other components so that when the full nonlinear equations were invoked the input 
vortex flow was dominated by the most unstable mode. A drawback of this method 
is that the lead-in time in which the linear mode develops proves to be a significant 
proportion of the total computation time. To alleviate this difficulty we used the 
program described by Bassom Otto (1992) to compute the unstable eigenfunction 
directly. This had the effect, of making the lead-in time to the nonlinear equations 
redundant and thus we could start our computations of the full nonlinear system almost 
immediately. Substantial reduction in the computational time was also achieved by 
using a stretched grid in the y-direction (normal to the cylinder surface) as opposed 
to a uniform one, and the scheme eventually chosen for this process was taken from 
Macaraeg, Streett Hussaini (1988). A grid is required which encompasses the region 
between y p = 0 and y p = y mar , where the subscript p denotes the physical coordinate 
and y ma x denotes some outer bound at which the asymptotic forms of the solution of 
the system (2.7) are supposed to have been attained. A notional computational grid 
0 < y c < 1 is introduced which is related to y p via 

Umax (f-2 1) Vc 

yp — / o \ c j i 

( c 2-y?) 

with the value (02 — 1) taking values between 0.2 and 1. The quantity (C 2 — 1) rep- 
resents the degree of stretching between the large and small steps whilst c\ controls 
the rate of stretching, 1 < C\ < 6. For most of the calculations reported here we used 
ci = 2.4 and C2 = 2 which allows 100 grid points to be distributed between y p = 0 
and y p = 50. The resulting distribution of points and the corresponding step lengths 
are illustrated in figure 1. Notice that as y — *■ y max the step lengths become greater 
than unity which in a finite difference scheme may induce errors. In the present case 
gradients of the dependent functions in this regime are minuscule and so this possible 
difficulty did not arise. This particular non-uniform stretching of the grid points al- 
lowed a twenty-fold reduction in the number of grid points over the regular grid used 
in Bassom & Hall (1991) and a four-fold improvement over the piecewise constant step 
length grid used in Bassom Otto (1992). 

As previously mentioned the solution strategy used here is essentially that de- 
scribed in Hall (1988). Suppose that the solution (U 0 , V 0 , U 0z , <f> 0 ) is known at some 
specified station, x, and suppose further that a guess is made for the solution at x + e. 
^£.q°\ V^°\ 0i O> )- The nonlinear terms on the right-hand sides of equations (2.7 a,b) 
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were calculated 

/ r (l) T-d) J\) 

y- 1 ' ' 1 • 

and the process repeated. W lien the difference between successive iterates was than 
some norm convergence was deemed to have occurred and this process marched onto 
the next step. For a computational grid composed of M points and a calculation with 
A modes retained then the convergence norm used was 


using these guesses and the system solved to provide updated values 
)• These updated variations were then put into the nonlinear terms 


fc=.V/2 — 1 j=M 

E E[iC + " 

k=-\ 7*2 >=1 


r r ( n ) I I T A 0 T 11 ) I 

L j.h I + i V j,k ~ V j,k I 


< 10~ 6 X. 


where £ Jj^ denotes the n th iterate of the A: th Fourier component of the flow quantity 
U evaluated at the j ih y- position. 

The equations (2.7 a,b) may be discretised for any particular mode into the forms 


i + 2 “I” brn C? 


+ dji 


•1 + 


+ /» 


- V n 


and 


(m = 3, * • * \I — 1 ) 
(3.1a) 


0 m ^ in + 1 ^ T7i ^ jn + l m U m -\ T ■ (m — 2, * * ’A/ — 1) (3.1&) 


This system was solved using a technique outlined in Appendix B of Bassom Otto 
(1992). This solves the problem (3.1) in one sweep rather than treating the pair (3.1 
a,b) as distinct pentadi agonal and tridiagonal problems. In essence, the equations are 
solved by performing an outward forward elimination followed by a back substitution. 

Earlier we alluded to the fact that significant computational saving was obtained 
by utilising Fast Fourier Transforms in preference to decomposing the flow quantities 
into their respective components. Standard procedures were implemented so that to 
change between physical and transform spaces we write 


K 2 




= X j€[l,JV], 


k=-Z- 
K 2 


1 j=N 

J'k = y X lkZj ' 



N 

2 


i], 


where Zj — 2i r(j — 1)/N for j G [l,iV] and where quantities with tildes are in the 
transform space and those without in the physical space. Denier k Hall (1992) were 
able to restrict themselves to using a cosine basis due to the nature of Gortler vortices 
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with two-dimensional boundary layers but the addition of temporal periodicity and 
crossflow prevented us from doing this. 

For the majority of the calculations presented below the number of modes retained 
was 16 and a fairly large step used in the streamwise direction, typically e = 0 . 01 . 
Denier Sc Hall (1992) found in their calculations that these parameter choices gave 
results to within graphical accuracy. Further testing, using 100 points in the normal 
direction and with y mar = 50 has confirmed that similar choices are satisfactory for 
the three-dimensional computations performed here. 


§4 Results 

We have detailed the numerical method by which we investigated the solution 
properties of system (2.7) although as yet we have left the definition of the amplitude 
of the initial vortex unspecified. As described, we initiated our computations with a 
multiple of the eigenfunctions of the linearised versions of system (2.7). For a specified 
vortex wavenumber k , frequency D and crossflow A the method outlined by Bassom 
Sc Otto (1992) was used to compute the corresponding linearised growth rates 3 r and 
the respective eigenfunctions normalised so that the energy defined by 



v 2 )dy 


is equal to A 2 . We refer to A as the amplitude of the initial condition. 

Our first calculations were performed primarily as a verification of our code against 
the established results of Denier Sc Hall (1992). Thence we considered the case of a 
purely two-dimensional boundary layer and steady vortices (A = F! = 0). Denier Sc 
Hall (1992) found that as the vortex evolves nonlinearly downstream then at some 
point their computations broke down in a singularity. Careful investigation revealed 
that at this location the skin friction was of the flow vanishes and then any solution 
scheme which relies upon a marching technique becomes invalidated. The vanishing of 
the skin friction was interpreted as being indicative of the vortices breaking away from 
the wall and moving into the core of the boundary layer. Our results for various initial 
amplitudes A are indicated on figure 2 where we show the location of the breakdown 
point Xb as a function of A for the most unstable Gortler mode with vortex wavenumber 
k — 0.476 (the disturbance was introduced at x = 1 so the distance travelled by the 
perturbation before breakdown is — 1). Not surprisingly, the breakdown location 
Xb is a monotone decreasing function of the amplitude A and as A — > 0 so Xb — > oc 
and the linear problem is retrieved. For all the other calculations reported upon here 
a similar trend is observed so that in all cases we chose A = 0.2; this selection was 
made purely for illustrative purposes and has no special significance whatsoever. 
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In figure 3 we demonstrate another feature of the two-dimensional results. Figure 
3a indicates the dependence of the linear vortex growth rati- .j r as a function of the 
vortex wavenumber k for the' first two modes. As is now well documented by Denier 
>:t al. (1991). Denier &: Hall (1992) the most unstable mode has growth rate 3 r = 
0.312 at k = 0.476 and ;i r — » 0 as k — > 0 or k -> oc. We note however, that for a 
significant range of scaled wavenumbers, roughly 0.3 < k < O.S. the vortex growth 
rate is not, much reduced from the maximum in as much as it is at least 00v< of 
the maximum. This will have important consequences for later results. Figure 3b 
illustrates breakdown point Xh as a function of k for the first two modes. As is to be 
expected for fixed initial amplitude A and fixed wavenumber the first mode always 
breaks away from the wall before the second one and this trend continues for higher 
modes. However, it is also observed that for the primary mode. x h is a monotone 
decreasing function of k. Therefore for any fixed initial amplitude it is not the most 
unstable mode which breaks up first. Indeed we previously remarked that for a whole 
range of wavenumbers surrounding the most unstable value linear growth rates are not 
too different from the largest growth rate. This has potentially important, consequences 
for a number of practical flows as it demonstrates that the breakdown properties of the 
flow are crucially dependent on the nature of the physical characteristics of the flow 
and are sensitive to the nature of the evolution of the flow. More precisely, suppose 
that the vortex motion starts with extremely small amplitude. Then one would expect 
there to be a large distance over which the motion develops essentially in a linear 
manner. During this time the most unstable linear mode would overwhelm modes of 
other wavelengths and once the vortices had grown sufficiently so that nonlinearity is 
important the flow behaviour would be dominated by that of the most unstable mode. 
On the other hand if the initial vortex amplitude is not tiny, nonlinear effect, are likely 
to be important a relatively short distance downstream by which stage it is unlikely 
that the most unstable mode dominates the others. In this case the behaviour of the 
most unstable vortex is not likely to dictate the properties of the breakdown of the 
flow. 

We turn now to consider cases with non-zero crossflows. As noted by Bassom 
Hall (1991) we may restrict our attention to cases in which the crossflow parameter 
A > 0 since by suitably transforming the system (2.7) we can relate flows with A < 0 
to appropriate counterparts with A > 0. In figure 4 we recall the results of increasing 
crossflow on the linear, stationary vortex mode. As discussed in detail by Bassom 
Hall (1991) the effect of crossflow on linear vortex structures of wavelength 0(G~s ) is 
primarily a stabilising one. Indeed figure 4 illustrates that when A > 0.410 the vortex 
mode is stabilised for all wavenumbers in the 0{G * ) regime. 


11 



Figures 5 a&; b show the variations of the maximum growth rate ($r) maz and 
the corresponding vortex wavenumber k maz for increasing crossflow parameter A. As 
already noted, increasing A lowers (3 r ) maz and it is the mode with wavenumber k ~ 
0.407 which is the last to be stabilised. Further k maz is not a monotone function of A 
as might have been anticipated. Figure 5b, as well as showing (d r ) maz , indicates the 
growth rate of the mode with wavenumber k = 0.476 (the value of k maz for the case of 
zero crossflow). It can be seen that the difference in growth rates of the most unstable 
mode and that with k — 0.4(6 is very small over the whole range of A considered. 
This reinforces the earlier comment that there is a significant wavenumber regime over 
which the linear growth rate of the vortex mode is almost constant. 

For each crossflow A we integrated the appropriate most unstable linear eigenmode 
with initial amplitude A = 0.2 until breakdown occurred, and then repeated the 
experiment with vortex wavenumber k - 0.476. The results are summarised in figure 
5c. It is seen that as well as exercising a stabilising influence on the linear mode, 
increasing the crossflow tends to have a similar effect on the nonlinear modes. Clearly 
the greater A, the greater the delay downstream before the initially most unstable 
mode breaks away from the wall. Additionally, for the same initial amplitude A and 
crossflow A the mode with k = 0.476 breaks up before the linearly most unstable mode. 
Further investigations have suggested that for all choices of A and vortex frequency 
ki then if two linear modes of different modes ki > k 2 are such that if initially have 
the same amplitude then if they are marched downstream then it is the one with the 
higher wavenumber which breaks down first. 

Finally, we discuss in more detail the influence of unsteadiness on our findings. In 
their study Bassom Sz Hall (1991) made some comments concerning the properties of 
time— dependent linearised vortices. For the majority of their work these authors were 
primarily interested in examining neutrally stable modes although they did compute 
a few non-neutral ones (see their figure 16). Bassom & Otto (1992) showed that for 
fixed wavenumber k then as the frequency Cl of the mode increases so the crossflow 
required to maintain neutral stability grows. In particular it was shown that for a 
non-dimensional vortex frequency Cl the stability properties of the vortex are sensitive 
to the sign of Cl. In figures 6 a& b we illustrate a facet of this sensitivity. For each 
frequency Cl and crossflow parameter A we show the wavenumber of the most unstable 
linear mode, k maz , together with the growth rate of that mode ( 0r) max ■ For Cl < 0, it 
it is observed that as A increases from zero so k maz decreases whereas for 0. > 0 this 
is not the case. Correspondingly, when Cl < 0 the growth rate of the most unstable 
mode is monotone decreasing in A whereas when Cl > 0 then as A increases from zero so 
there is a crossflow range over which (0 r ) maz increases. This increase is not indefinite 
however and there is a critical A, dependent upon Cl, after which the growth rate 
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decreases. We also note that if 0\ > O 2 > 0 then max^ ft r Aj < max^ 3 r A j 
so that over all Cl and all A > 0 the mode with the greatest growth rate is stationary 
and exists in a two-dimensional boundary layer. 

The breakdown characteristics of unsteady flows are described by figure 6c. For 
each frequency Cl and A we marched the linearly most unstable mode of amplitude 
A = 0.2 from x = 1. We can observe the somewhat conflicting roles played by 
crossflow and frequency. 

In the main, for a prescribed Cl increasing A delays breakdown whereas for fixed 
A and increasing Cl this phenomenon is enhanced. Notice, however, one important 
feature which runs against this general rule. For positive Cl then a small to moderate 
crossflow actually tends to promote breakdown although larger crossflows do reverse 
this effect. An attempt was made to verify this trend by considering larger values of 
Cl than those illustrated in figure 6. However problems were encountered as Cl grew 
and these difficulties can be attributed to a number of causes. Following on from work 
elucidated in Bassom k Hall (1991) it is the case that for small crossflow the most 
unstable linear mode first has a small wavenumber relative to the implied scaling. As A 
grows then the most unstable mode corresponds to an eigenfunction that moves away 
from the wall at y = 0. At Cl greater than about 2 this movement occurs quickly for 
small changes in A so that for quite moderate values of A the eiegnmode is far removed 
from the boundary. As found both by Bassom k Hall (1991) and Bassom k Otto 
(1992) the numerical solution of the governing equations becomes noil-trivial as the 
vortex moves out since boundary conditions need to be imposed at the wall. Charlv 
for modes concentrated away from the wall large changes in 3 , A or Cl can lead to 
almost imperceptible changes in the values of the eigenfunctions at the wall and thus 
reliable numerical convergence is rendered very difficult. However, our limited further 
computations for Cl > 1.5 are in agreement with the general behaviour described above. 

§5 Conclusions k Discussion 

In this w r ork we have detailed the nonlinear spatial evolution of unstable Gortler 
modes in a three-dimensional boundary layer. In particular, the roles played by vortex 
wavenumber, frequency, and the crossflow component of the underlying base flow have 
been described. We feel that of particular importance is our finding that ( all other 
factors being equal) of two modes of wavenumber within the 0{G*) regime the one 
with the smaller wavelength will be the first to breakdown. This then suggests that in 
practical situations it may not be the most unstable linear mode which is of ultimate 
importance. 

In many cases the relevant calculation to describe the breakdown of a flow is one 
of a receptivity type. In this scenario small disturbances within the boundary layer 
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or on the wall of the cylinder can trigger Gortler vortices and the precise method 
of this triggering frequently excites modes of a preferred wavenumber. If this occurs 
in practice then our calculations provide a description of the evolution of the mod*'. 
Conversely, if a range of wavenumber modes is excited two eventualities would seem to 
be possible. First, suppose that the initial disturbance is very small. Then it is to be 
expected that the perturbation travels a long way downstream before nonlinearity has 
significant effect, the most unstable linear mode will be dominant before this point and 
its evolution characteristics will essentially describe that of the whole flow. Second, 
suppose that the initial perturbation is not so small. Our results summarised bv figure 
3 have shown that although there is a unique most unstable mode for each crossflow A 
and frequency Cl vortices with wavenumbers in a fairly-large region surrounding that 
of the most unstable mode have growth rates not very different from the maximum. 
Therefore, by the time nonlinearity is significant it is not clear that the most unstable 
mode would be dominant and the breakdown characteristics of the whole flow would 
involve calculations more involved than those reported here. However we have shown 
that for given mode amplitude it is the higher-wavenumber modes which appear to 
breakdown first so that these components of a spectrum of excited modes may well 
prove to be the important ones. 

Denier & Hall (1992) showed that when their calculations for nonlinear modes in 
two-dimensional boundary layers terminated, this corresponds to the skin friction of 
the flow vanishing at some point. Once this happens marching schemes as used both 
here and in Denier & Hall (1992) cannot be continued. We confirm this finding for 
our three-dimensional cases as well but we also observed that before the skin friction 
\ anishes the velocity profile develop inflection points at positions awav from the wall. 
The appearance of those inflection points suggest that the flow will become susceptible 
to Rayleigh waves which would give an alternative route to the ultimate breakdown. 
The analysis of these modes would be of interest. 

Finally, we recall that all our calculations have been concerned with considering 
the evolution of perturbations of a specified wavenumber. Of course in some situations 
a spectrum of modes may well be present. W r e have identified situations in which we 
might expect one mode to dominate the others before nonlinearity sets in but in the 
other cases calculations would be needed which account for an initial perturbation 
which contains a number of modes. The development of a code to perform such calcu- 
lations might well be formidable but it would give the definitive theoretical description 
nonlinear Gortler vortex behaviour in three-dimensional boundarv layers. 
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Figure 1: Distribution of grid points and corresponding spacing for the calculations 
performed here. The computational coordinate y r satisfies 0 < y c < 1 ail( l related 
to the physical coordinate y p by y p — 50y c /(2 — y z c ) 1A ■ 
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Figure 3: 

a) Linear vortex growth rates 8 r as a function of wavenumber k for the two most 
dangerous modes in a two-dimensional boundary layer. 

b) Breakdown point x b as a function of wavenumber k for the modes in figure 3a with 
assumed initial amplitude A = 0.2. The symbols denote corresponding points on the 
curves of a) and b). 
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Figure 5: 

a) Wavenumber of the most unstable stationary mode and b) its corresponding growth 
rate as a function of crossflow A. 

c) Breakdown point as a function of crossflow A for z) the most unstable mode and 
ii) the mode with k = 0.476 for initial amplitude A = 0.2. 
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Figure 6: 

a) Wavenumber of the most unstable nonstationary modes and b) corresponding 
growth rates f3 r as functions of A. Here we have considered Cl = —1, 0.5, 0, 0.5, 1, 1.5. 
c) Breakdown point xj as a function of crossflow A for the modes with frequencies Cl 
and wavenumbers k as given in figure 6a) and initial amplitude A = 0.2. 





